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By analyzing the temperature (T) and density (n) dependence of the measured conductivity (a) 
of 2D electrons in the low density (~ 10 11 cm -2 ) and temperature (0.02 - 10 K) regime of high- 
mobility (1.0 and 1.5 x 10 4 cm 2 /Vs) Si MOSFETs, we establish that the putative 2D metal-insulator 
transition is a density-inhomogeneity driven percolation transition where the density-dependent 
conductivity vanishes as er(n) oc (n — n p ) v , with the exponent p ~ 1.2 being consistent with a 
percolation transition. The 'metallic' behavior of o~{T) for n > n v is shown to be well-described by 
a semi-classical Boltzmann theory, and we observe the standard weak localization-induced negative 
magnetoresistance behavior, as expected in a normal Fermi liquid, in the metallic phase. 
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The so-called two-dimensional (2D) metal-insulator 
transition (MIT) has been a subject 1 -! 2 - of intense activity 
and considerable controversy ever since the pioneering 
experimental discovery 3 of the 2D MIT phenomenon in 
Si MOSFETs by Kravchenko and Pudalov some fifteen 
years ago. The apparent MIT has now been observed 
in almost all existing 2D semiconductor structures, in- 
cluding Si MOSFETs&i, electrons^ and holesS&i&i 1 . 
in GaAs/AlGaAs, and electrons in Si/SiG o 12 ' 13 systems. 
The basic phenomenon refers to the observation of a car- 
rier density-induced qualitative change in the tempera- 
ture dependence of the resistivity p{n,T), where n c is 
a critical density separating an effective 'metallic' phase 
(n > n c ) from an 'insulating' phase (n < n c ), exhibiting 
dp/dT > (< 0) behavior typical of a metal (insulator). 

The high-density metallic behavior (n > n c ) of- 
ten manifests in a large (by 25% for electrons in 
GaAs/AlGaAs heterostructures to factors of 2-3 in Si 
MOSFETs) increase in resistivity with increasing tem- 
perature in the low temperature (0.05 - 5 K) regime 
where phonons should not play much of a role in resis- 
tive scattering. The insulating regime, at least for very 
low (n <C n c ) densities and temperatures, seems to be 
the conventional activated transport regime of a strongly 
localized system. The 2D MIT phenomenon occurs in rel- 
atively high-mobility systems although the mobility val- 
ues range from 10 4 cm 2 /Vs (Si MOSFET) to 10 7 cm 2 /Vs 
(GaAs/AlGaAs) depending on the 2D system under con- 
sideration. The 2D MIT phenomenon is also considered 
to be a low-density phenomenon although, depending 
on the 2D system under consideration, the critical den- 
sity n c differs by two orders of magnitude (n c ~ 10 11 
cm~ 2 in 2D Si and ~ 10 9 cm -2 in the high-mobility 
GaAs/AlGaAs heterostructures. The universal features 
of the 2D MIT phenomenon are: (1) the existence of a 
critical density n c distinguishing an effective high-density 
metallic {dp/dT > for n > n c ) phase from an effective 
low-density insulating {dp/dT < for n < n c ) phase; and 
(2) while the insulating phase for n < n c seems mostly to 



manifest the conventional activated transport behavior, 
the metallic temperature dependence at low {n > n c ) 
densities is universal in the sense that it manifests a 
very strong temperature dependence, not seen in stan- 
dard 3D metals where p{T) is temperature independent 
in the T < 10 K Bloch-Griineisen transport regime. 

The excitement and controversy in the subject arise 
from the deep conceptual questions associated with the 
nature of the MIT and the metallic phase. In particular, 
it is theoretically well-established 1 ' that a non-interacting 
(or weakly interacting) disordered 2D electron system is 
an insulator at T — 0, and therefore it follows 1 that if the 
2D MIT is a true T — quantum phase transition (QPT), 
then the 2D metallic phase, if it survives all the way to 
T = 0, must necessarily be a novel and exotic non- Fermi 
liquid phase since it cannot be connected adiabatically 
to the non-interacting 2D electron system, which is al- 
ways localized in the presence of any disorder. The alter- 
native possibility 2 " is that the metallic phase is allowed 
only at finite temperatures, and the 2D MIT is not a 
QPT, but a density-induced crossover from a weakly lo- 
calized 'effective' metallic phase to a strongly localized 
insulator. These two qualitatively distinct viewpoints 
have both received support in the literature, and experi- 
mental results have been presented claiming support for 
the QPT™ and the finite-temperature crossoveriii view- 
points, respectively Whether the observed 2D MIT is a 
novel (T = 0) QPT leading to an exotic 2D non-Fermi 
liquid metallic phase established by interaction effects 
or is a finite-temperature crossover between an effective 
metallic phase and the insulating phase is obviously a 
question of great fundamental importance. 

Previous measurements of 2D MIT phenomena in Si 
MOSFETs have been discussed in terms of quantum crit- 
ical phenomena 1 *!^. However, for 2D electrons and holes 
in high mobility GaAs/AlGaAs heterostructures, it has 
been shown that the 2D MIT can be described as a per- 
colation transitio n 7 ' 11 ! 15 . Thus, a natural question to ask 
is whether the 2D MIT in Si MOSFETs is really different 
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than for GaAs systems, or whether both Si and GaAs 2D 
systems undergo the same universal percolation-driven 
transition. The Si MOSFET is different than GaAs sys- 
tems in that it has much lower mobility and the disorder 
in Si MOSFETs is thought to be dominated by 'short- 
range' scattering, which is potentially different than for 
other 2D systems, such as in GaAs where 'long-range' 
scattering may play a larger role. This point has been 
recently emphasized in the literature and it has been pro- 
posed that the physics behind the 2D MIT for Si MOS- 
FETs versus GaAs systems may differ due to differences 
in mobility and range of the disorder potential (short vs. 
long-range scattering)ii. 

In this article, we provide the first experimental ev- 
idence that the MIT for Si MOSFETs is a finite- 
temperature, density inhomogeneity driven semi-classical 
percolation transition, and is therefore not a QPT. Our 
work shows that the 2D MIT is a universal percolation 
transition for both GaAs and Si, in contrast to claims 
of a QPT in SiH As mentioned earlier, the 2D MIT 
in high-mobility GaAs electron and hole systems has al- 
ready been experimentally demonstrated to be a perco- 
lation transition ^ 11 ' 15 , but our work provides the first 
compelling evidence supporting the percolation scenario 
in Si MOSFETs. We also establish that the metallic 
phase is a conventional Fermi liquid by showing that (1) 
the measured transport properties in the effective metal- 
lic phase qualitatively agree very well with the conclu- 
sions of a semi-classical Boltzmann transport theory^ 
taking into account the detailed density and tempera- 
ture dependence of carrier screening properties, and (2) 
the expected negative magnetoresistance signature of the 
conventional weak localization behavior— is clearly man- 
ifest experimentally in the effective metallic phase of our 
samples. 

The two samples used in this study are Si MOSFET 
structures with a peak mobility of ~ 1.5 x 10 4 cm 2 /Vs 
(sample A) and ~ 1 x 10 4 cm 2 /Vs (sample B). For both 
samples, the 2D electron gas (2DEG) resides at the Si- 
Si02 gate oxide interface, where the Si02 thickness is 
nominally 35 nm for sample A and 10 nm for sample B. 
The gate oxide is thermally grown in dry O2 at 900° C on 
FZ (lOO)-oriented high-resistivity (> 5 kf2-cm for sample 
A and > 100 17-cm for sample B, at room temperature) 
p-type Si wafers. Ohmic contacts to the 2DEG consist of 
n + Si regions formed by implantation of As; these con- 
tacts are not self-aligned, but are instead patterned and 
implanted before the gate oxide growth. An n + polySi 
gate, patterned to form a 15 /im x 100 /im Hall bar, is 
used to induce carriers. The 2DEG resistivity is exper- 
imentally determined via standard four-terminal lock-in 
measurements and the density is calibrated via measure- 
ments of the Hall resistivity. 

In Fig. 1 we show the density-dependent mobility of 
samples A and B along with the theoretically calculated 
mobility curves. Our theoretical calculations 1 ? include 
the temperature-dependent finite wavevector screening of 
the charged impurity Coulomb scattering and the short- 



2 




(a) " 


Q 1 .5 

CM 


/ \ 


- 


1 1-0 






E0.5 


I / B 










0.0 







im 


(t)) ■ 




• ll SR 





012345012345 



nrjl -2\ 

n (10 cm ) 



n (10 12 cm~ 2 ) 



FIG. 1: (color online) (a) Experimental mobility /i as a func- 
tion of density n for two samples, A and B, at a temperature 
T = 0.25 K. (b) Theoretical mobility as a function of density. 
High (low) density mobility is limited by the surface roughness 
Msr (charged impurity fj,i m ) scattering. Different curves cor- 
respond to different charged impurity densities (rn = 6, 9, 12, 
and 24 x 10 10 cm -2 from top to bottom), but with the same 
surface roug hness (A = 1.9 A, A = 26 A). 



range surface roughness scattering at the Si-Si0 2 inter- 
face. It is well-known that the long-range charged im- 
purity scattering and the short-range surface roughness 
scattering dominate respectively the low and the high 
carrier density regimes of transport in Si MOSFETs 19 . 
We show Fig. la and b together to allow the reader 
to judge the level of agreement between experiment and 
theory, which depends on several scattering parameters 
(e.g. the density and spatial distribution of charged im- 
purity centers in the sample and the height and lat- 
eral correlation functions of the interface roughness func- 
tion). The mobility (fj,) first rises steeply with density 
(n) reaching a maximum at a characteristic density n m , 
beyond which /z falls off slowly with increasing n, due 
to short-range surface roughness scattering at the Si- 
Si02 interface^. At low carrier densities (n < 10 12 
cm~ 2 ) however, long-range Coulomb scattering by unin- 
tentional random charged impurities invariably present 
on the insulating oxide side of the Si-Si02 interface dom- 
inates the 2D mobility, and, as has been emphasized 20 
in the literature, all Si 2D MIT behavior manifests it- 
self at low carrier densities (n ~ 10 11 cm -2 ) where the 
charged impurity scattering dominates transport and the 
short-range interface roughness scattering is negligible. 
We emphasize the fact that the 2D MIT phenomenon 
(i.e. n c and the density range for the 'anomalous' metal- 
lic phase with a strongly temperature dependent 2D re- 
sistivity) occurs in the density regime n ~ 10 11 cm -2 
which, being substantially below n m , is completely dom- 
inated by the long-range screened Coulomb scattering. 
This dominance of Coulomb scattering in the 2D MIT 
physics is crucial for the percolation physics which arises 
from the nonlinear failure of screening of the long-range 
Coulombic charged impurity potential, to be discussed 
below. The short-range surface-roughness scattering, 
which dominates MOSFET transport for n > 10 12 cm" 2 , 
plays no role in the 2D MIT phenomena since both the 
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FIG. 2: (color online) (a) Experimental resistivity p in 
units of h/e 2 as a function of temperature for sample 
A at 2D electron densities (from top to bottom) n = 
1.07, 1.10, 1.13, 1.20, 1.26, 1.32, 1.38, 1.44, 1.50, 1.56, 1.62, and 
1.68 x 10 11 cm -2 . Inset (b) Enlarged view of data in the effec- 
tive metallic regime at n = 1.56, 1.62, and 1.68 x 10 11 cm -2 . 
Solid lines are guides to the eye. (Data below 0.1K may not 
be reliable due to electron heating.) (c) Theoretically calcu- 
lated temperature and density dependent resistivity for sam- 
ple A for densities n = 1.26, 1.32, 1.38, 1.44, 1.50, 1.56, 1.62, 
and 1.68 x 10 11 cm -2 (from top to bottom). Inset shows the 
resistivities for n = 1.56, 1.62, and 1. 



critical density for the transition and the density regime 
where the strong metallic behavior (i.e. the strong tem- 
perature dependence of the 2D conductivity) is observed 
are well below the operational regime of the short-range 
surface roughness sacattering. 

In Fig. 2 we show our measured (Fig. 2(a) - (b)) 
and calculated (Fig. 2(c) - (d)) temperature and den- 
sity dependent resistivity (p(T, n)) for sample A with the 
maximum mobility of 1.5 x 10 4 cm 2 /Vs (see Fig. 1(a)). 
The results for sample B are similar, but with a larger 
value of the critical density, n c ~ 2.5 x 10 11 cm -2 (and 



somewhat weaker temperature dependence), consistent 
with its lower (and consequently, higher disorder) peak 
mobility ~ 1 x 10 4 cm 2 /Vs. The classic 2D MIT be- 
havior is apparent in Fig. 2(a) where p(T) for various 
densities manifests the clear distinction between metallic 
(dp/dT > 0) and insulating (dp/dT < 0) behavior sep- 
arated (visually) by a critical density n c ~ 1.4 x 10 11 
cm~ 2 . We note that the temperatures quoted in Fig. 2a 
and b are those of the cryostat cold finger; although the 
resistivity continues to evolve below T ~ 100 mK (albeit 
slowly), we cannot claim in this regime (T < 100 mK) 
that the 2D electrons and the cold finger are in thermal 
equilibrium down to the base temperature of our cryo- 
stat (20 mK). The inset in Fig. 2(b) shows in detail 
the metallic behavior for n > n c , which manifests more 
than a factor of 2 increase in p(T) for T ~ 0.1 — 7 K 
before decreasing slightly at higher temperatures. The 
strong initial increase in p(T) with T arises^ from the 
temperature-induced weakening of screening at these low 
densities as T/Tp increases from 0.01 at low tempera- 
tures essentially to 1 at T ~ 8 K since Tp « 7.3 K for 
n = 10 11 cm~ 2 . (Phonon scattering is unimportant in Si 
MOSFET for T < 10 K.) 

Finally, in Fig. 2(c) we show our theoretically calcu- 
lated p(T, n) where the Boltzmann transport is calcu- 
lated using the screened charged impurity scattering as 
the only resistive scattering mechanism 1 ^. Although the 
ratio of the Coulomb to Fermi energy, r s , can be as large 
as ~ 15 near the MIT for this sample, suggesting that 
Coulomb interactions could play an important role, we 
note that the basic features of the experimental metal- 
lic phase are well-captured by the screening theory with 
p(T) rising approximately by a factor of 2 initially and 
then decreasing slightly at higher (lower) temperatures 
(densities) around T /Tp > 1 due to quantum-classical 
crossover—. The theoretical results, following the work 
of Refji&, are being shown here in order to demonstrate 
that a realistic semiclassical transport theory including 
the screened long-range scattering by charged impurities 
completely captures all the main features of the experi- 
mental results with the basic physics being a strong tem- 
perature and density dependent modification in the ef- 
fective screened disorder due to the large change in the 
effective T/Tp in the experimental T ~ 0.1 — 10 K tem- 
perature range. We point out that, by definition, the 
theoretical transport results in Figs. 2(c) and (d) ap- 
ply only in the metallic regime since the metal-insulator 
transition cannot be captured in a Boltzmann theory. 

Having established the phenomenology of the 2D MIT 
behavior in Fig. 2, we now come to the nature of the den- 
sity driven transition itself. We plot in Fig. 3 our mea- 
sured conductivity cr(n) as a function of n to see if a per- 
colation behavior, cr(n) ~ (n — n p ) p , where n p , p are the 
percolation transition density and exponent respectively, 
is manifested. Our fit yields p ~ 1.20, n p 1.2 x 10 11 
cm~ 2 for sample A and p ~ 1.24, n p sa 2.0 x 10 11 cm" 2 
for sample B. The fits describe cr(n) well over nearly two 
orders of magnitude change in conductivity. The devia- 
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FIG. 3: (color online) (a) Main plot: Points show experimen- 
tal conductivity a in units of e 2 /h versus electron density n 
for sample A at T = 0.4 K (closed symbols) and B at T = 0.25 
K (open symbols). The solid lines are fits to the data of the 
form <r = A(n — n p ) p . The upper and lower insets show the 
exponent p and critical density n p , respectively, as a function 
of temperature for sample A. Solid lines are a guide to the 
eye. The open symbol in the upper inset shows p for sample 
B at T = 0.25 K. 
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FIG. 4: (color online) Main plot: Points show n p vs. 1/T 
obtained from a percolation fit to the a vs. n data performed 
by holding the exponent p — 1.31 fixed (see text). The solid 
trace is a fit of the form n p = no + Ae~ b ^ T . Inset: Points 
show experimental data a vs. n for sample A at T = 0.4 K, 
while the solid line is a percolation fit with fixed p = 1.31. 



tions of the fit from the data seen at the lowest densities 
( n ^ n c) are expected since our percolation analysis ap- 
plies only in the metallic regime. At high densities the fit 
is also expected to deviate since the percolation fit is a 
critical scaling behavior which is expected to be best for 
densities close to n p , and also the role of surface rough- 
ness scattering becomes increasingly important at high 
densities. The fact that our best fit values of n p and 
p for sample A are almost independent of temperature 
at low temperatures is a good indicator for the percola- 
tion transition being the correct description for the 2D 
MIT phenomenon shown in Fig. 2. Our numerically ex- 
tracted best- fit exponent p w 1.2 is sufficiently close to 
the percolation conductivity exponent of 1.31 for us to 
feel confident that the observed 2D MIT in Si MOSFETs 
is a density inhomogeneity-driven percolation transition. 

In Fig. 4, we show an alternate (and perhaps more ap- 
propriate, if the transition is indeed a percolation tran- 
sition) method of extracting the critical density from a 
percolation fit, where we hold both the exponent p and 
the prefactor A fixed, and allow only n p to vary. We 
use p = 1.31, the value expected for percolation in 2D, 
and determine A from the percolation fit to the T = 0.4 
K data, as shown in the inset to Fig. 4, also holding 
p = 1.31 fixed. If the 2D MIT is indeed a percolation 
transition, then this method of obtaining n p is more ap- 
propriate than allowing p to vary since the exponent is 
universal and the critical density is not. As seen in the 
inset to Fig. 4, the fit to our a vs. n data with p = 1.31 
fixed is very good. As can be seen from the main figure, 
the extracted n v obtained via this method decreases with 



decreasing temperature, eventually saturating at the low- 
est temperatures. A fit of the form n p = uq + Ae~ h l T is 
shown as the solid line in the main plot of Fig. 4. The fit 
yields n p = 1.04 + 1.62e~ 4 20/T . The parameter b = 4.2 
K is an effective "energy gap" that can be calculated if 
the impurity distribution is known. We note that this 
implies a T = MIT transition density of 1.04 x 10 11 
cm , which is considerably less than the nominal crit- 
ical density of 1.4 x 10 11 cm -2 one would have inferred 
from the temperature-dependent resistivity behavior in 
Fig. 2(a). The fact that the temperature-dependent re- 
sistivity by itself gives a strong over-estimate of the 2D 
MIT critical density is already known in the literature. 

In Fig. 5 we show our observed weak localization be- 
havior by plotting the expected^ negative magnetoresis- 
tance in the metallic phase as a function of an applied 
weak magnetic field. The weak localization data in Fig. 
5 is fitted to the standard di-gamma function behavior 
expected of a disordered 2D system^ 

-=-4?[*G+?)-*G+?)]' w 

where tb = h/AeB^D and D is the diffusion coefficient. 
The behavior is quite normal, explicitly demonstrating 
that our observed metallic phase (n > n Cl n p ) is indeed 
the usual weakly localized Fermi liquid metallic phase, 
and not some exotic non-Fermi liquid T = metal. The 
observation of the expected weak localiation behavior in 
the 2D metallic phase along with the percolative nature 
of the transition is strong evidence that the experimental 
2D MIT is not a QPT, but a crossover. 
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FIG. 5: (color online) Measured weak localization correction 
to the resistivity, Ap xx /pQ in units of e 2 /h (with offset for 
clarity) versus perpendicular magnetic field B at densities n = 
1.51, 1.45, and 1.63 x 10 11 cm' 2 at T = 100 mK. The red lines 
are fits to the standard weak localization correction (see text). 
The momentum relaxation times r are determined from the 
mobility, while the weak localization fits yield a correction 
amplitude a = 0.25 and phase relaxation time = 33, 32, 
and 28 ps (top to bottom). 

Before concluding, we briefly discuss the physics under- 
lying the percolation transition. As emphasized above, 
the 2D MIT phenomenon occurs in the transport regime 
dominated by long-range Coulombic charged impurity 
scattering. It was pointed out some years ago— that 
such a transport regime is susceptible at low carrier den- 
sities to a semiclassical percolation MIT since large scale 
density inhomogeneities (i.e. puddles) could appear in 
the system due to the nonlinear failure of screening of 
the charged impurities at low carrier densities. One can 
approximately estimate 21 the percolation transition den- 
sity, n p , by considering the inhomogeneous density fluc- 
tuations in the system, leading to n p ~ 0.2-y/nt/d, where 
rii is the effective 2D charged scatterer density and d 
is the effective separation between the 2D carriers and 
the scatterers. Using m = 6 and 11 x 10 10 cm" 2 , the 
value needed to get the correct mobility at higher den- 
sities for sample A and B, respectively (these values of 
rii are roughly consistent with C-V measurements which 
estimate the density of fixed charge near the Si-Si02 in- 
terface to be ~ 8 x 10 10 cm" 2 ), and taking d ~ 5 — 6 
nm, which is consistent with the 2D electron layer in 
the Si MOSFET being about 5 nm away from the in- 
terface due to the self-consistent potential produced by 
the electrons themselves, we get n p ~ 0.9 and 1.2 x 10 11 
cm" 2 for sample A and B, respectively, which is close to 
our experimentally extracted percolation density in Fig. 
3(a). (We note that d corresponds to the effective separa- 
tion between the 2D electrons and the charged impurities, 



which must incorporate the fact that the 2D electrons are 
on the average about 5 nm inside Si at n ~ 10 11 cm" 2 .) 
The 2D conductivity percolation exponent is numerically 
known to be around 1.31, which is close to the exponent 
value (1.20 for sample A, 1.24 for sample B) we get in our 
analysis. Finally, we note that the puddles produced by 
the nonlinear failure of screening of charged impurities 
in low-density 2D systems have been directly observed in 
2D G&AtM systems. 

In conclusion, we have experimentally established that 
the 2D MIT in Si MOSFETs is a percolation-induced 
transition driven by the dominance of the long-range 
charged impurity disorder. At low carrier densities, 
the failure of linear screening leads to the formation 
of inhomogeneous puddles in the 2D density landscape, 
which then produces a semi-classical percolation transi- 
tion. The nature of the percolation transition and the 
effective metallic phase is the same in Si MOSFETs and 
2D GaAs electron! and holeii systems. The metallic 
temperature dependence arises from the strong tempera- 
ture dependence of screening and the percolation transi- 
tion arises from the low-density failure of screening. The 
2D MIT is therefore a quantum crossover phenomenon, 
not a QPT, from an effective (weakly localized) metal- 
lic phase to a strongly localized insulating phase. The 
agreement between the percolative conductivity and our 
experimental work as well as our finding of the usual 
weak localization effect in our data conclusively rules out 
quantum criticality in the 2D MIT phenomena. Our work 
establishes that the 2D MIT is a universal percolation- 
driven transition, caused by the inhomogeneous density 
landscape induced by charged impurities at low carrier 
densities, in both GaAs-based systems ^ 11 ' 15 as well as 
in Si MOSFETs (our work). Screening of the charged 
impurity scattering produces the strong temperature de- 
pendence of the apparent metallic phase since low carrier 
densities automatically imply large values of T/Tp, and 
the eventual failure of screening at still lower densities, 
where there are just too few carriers to effectively screen 
the impurity charges, leads to the percolation MIT since 
the 2D system becomes inhomogeneous due to the for- 
mation of puddles. 
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